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Abstract. Building on results from previous studies of Cygnus X-1, we analyze Rossi X-ray Timing Explorer (RXTE) data 
taken when the source was in the soft and transitional spectral states. We look at the power spectrum in the 0.01 - 50 Hz range, 
using a model consisting of a cut-off power-law and two Lorentzian components. We are able to constrain the relation between 
the characteristic frequencies of the Lorentzian components, and show that it is consistent with a power-law relation having the 
same index (1.2) as previously reported for the hard state, but shifted by a factor ~ 2. Furthermore, it is shown that the change in 
the frequency relation seen during the transitions can be explained by invoking a shift of one Lorentzian component to a higher 
harmonic, and we explore the possible support for this interpretation in the other component parameters. With the improved 
soft state results we study the evolution of the fractional variance for each temporal component. This approach indicates that 
the two Lorentzian components are connected to each other, and unrelated to the power-law component in the power spectrum, 
pointing to at least two separate emission components. 
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1. Introduction 

Cygnus X-1 is one of the most studied X-ray sources, and is 
often quoted as the prototype black hole binary system. The 
source exhibits two main spectral states, commonly referred to 
as hard an d soft, with a brief intermediate state during tr ansi- 
tions (e.g., IZdziarski et alJ2002tlz'dzia rski & Gierlifiski^2004). 
Recently, a study of the broad-band spectra of Cyg X- 1 , mainly 
in the hard state, has been presented bv llbragimov et al.l (120051) . 
and a comprehensive study of all states of t he source in th e 
3-200 keV range has been carried out by Wil ms et al 1 (l2006l) . 
Several models have been proposed to explain the observed 
states and transitions. The two main components of such mod- 
els are usually a geometrically thin, optically thick accretion 
disk and a hot inner flow or corona. The models vary in the 
geometry and properties of mainly the corona/Comptonizing 
region (for a review on coronal models, especi ally in re- 
gard to timing characteristics, see IPoutanerilbOOll) . However, 
changes in the inner radius of the accretion disk are a com- 
mon source of variabili t y in many models (e.g ., Poutanen et dj 
1997'; 'E sin et al.lfl998l IChurazov et al J 1200 it [zdziarski et dl 



20021) . 

The changes in spectral state are evident also in the power 
density spectrum (PDS). In the hard state the PDS is charac- 
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terized by a flat-topped component up to ~ 0.2 Hz, where 
it breaks to a slop e that steepens at a few Hz (e.g., 
Belloni & Hasinger 199d lNowak et alJll999l) . In flie soft state 
the PDS is dominated by a f~ ' component that steepens around 
10 Hz (e.g.. lCui et alJl997ah . The 1996 state transition also re- 
vealed an intermediate PDS, showing a ~ / ' slope at lower 
frequencies, a flatter component around 0.3 - 3 Hz, and above 
that similar to the hard state PDS. In many observations there 
is evidence for a qtiasi-periodic oscillatioii (QP O) at a few Hz 
(iBelloni et al.ll996HCui et alJl997htlCui Il 999l). 

Previous studies of the PDS of Cyg X-1 have shown cor- 
relations between different temporal features as well as be- 
tween temporal and spectral components (e.g.. iGilfanov et aP 
Il999l) . Correlations between features in the hard st a te PDS 
were first reported in lWiinands & vari derKlisl(ll999l) . lNowakl 
fcoOOJ showed that the hard state PDS of Cyg X-1 could 
be well fit using Lorentzian components. This was expanded 
in a l arge scale study of the hard state PDS conducted by 
IPottschmidt et al. t2003.) . who also included some 'failed state 
transition' PDS. A systematic study of the PDS for al l spec- 
tral states of Cyg X-1 was conducted bv lAxelsson et alJ (l2005[ 
hereafter Paper I) using archival data from the Rossi X-Ray 
Timing Explorer (RXTE) satellite. Using a model consisting of 
a cut-off power-law and two Lorentzian profiles, we were able 
to fit the majority of PDS, and follow the components from 
hard state through the transitions and back. While the results 
showed that the Lorentzian components were present in a sig- 
nificant fraction of the soft state PDS, their behavior was not 
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well constrained. In this paper we reanalyze the soft state data 
and expand it with observations from an additional soft state 
in May- August of 2003, recently made available. With the im- 
proved soft state results we are able to decompose the PDS and 
study the contribution from each component separately. 

We begin this paper with a brief description of the 
data used, the analysis procedure and the model components 
(Sect. This is followed by a presentation of our results in 
Sect. |3l In Sect. |3 we discuss our findings, and put them in 
context of other results. Finally our findings are summarized in 
Sect.ia 

2. Observations and Data Analysis 

The data and method s of ana lysis used in this paper are in 
essence the same as in iPaoer j Therefore we only give a brief 
summary of the data and extraction process, with emphasis on 
the new data and the parts of the analysis that have been im- 
proved. For more details we refer to iPaoer It 

2.1. The archival data 

Since its launch on December 30, 1995, the RXTE satellite has 
observed Cyg X-1 in more than 750 pointed observations, with 
a typical on-source time of 3 - 4 ks. During this time, the All- 
Sky Monitor (ASM) instrument on board the satellite has also 
provided a continual coverage of the source in the 2-12 keV 
range. 

As in our previous study, the d ata are obtained w ith 
the Proportional Counter AiTay (PCA, Ijahoda etal ?1996') on 
RXTE, which consists of five identical Proportional Counter 
Units (PCUs), and is sensitive in the 2-60 keV energy range. 
The new soft state roughly spans the time from May 2003 to 
August 2003. We include all pointed observations made of the 
source during these four months, ~ 4 days per month. Note tha t 
some of these data are presented in Pottschmi dt et all (|2005). 
The data mode used was the 'Generic Binned' mode. We also 
calculate count (hardness) ratios, using data from the Standard2 
configuration. 

Lightcurves were extracted with the standard RXTE data 
analysis software FTOOLS, version 5.2, using standard screen- 
ing criteria: a source elevation > 10°, a pointing offset < 0?01 
and a South Atlantic Anomaly exclusion time of 30 minutes. 
The energy range of the data used for the PDS is ~ 2 - 9 keV, 
with the soft and hard bands of the count ratio taken as 2-4 keV 
and 9-20 keV respectively. For the PDS, the lightcurves were 
extracted with a time resolution of 10 ms. The lightcurves used 
in determining count ratios were extracted with 16 s time res- 
olution. All lightcurves were normalized to one average PCU 
using the FTOOL correctlc. 

2.2. Calculating the power spectra 

The modelling of the soft state PDS in ' Paper j (e.g.. Fig. 9) 
showed the need to extend the studied frequencies above 
25 Hz. This requires proper treatment of dead-time effects, and 
we therefore correct the estimated Poisson level. We u se the 
correction for general paralyzable dead-time of .Zhan^ et alJ 



(I1995') as presented in the first coiTective term in Eq. 6 of 
Jerniaan et al. (2000). We chose an upper limit of 50 Hz. 

To calculate the power density spectra the high-resolution 
light curves were divided into segments of 2'^* bins (~ 82 s). 
The resulting PDS were then gathered into groups of 10 and 
each group averaged into one final PDS. Each PDS thus cov- 
ers ~ 15 minutes of observation time, allowing changes oc- 
curing on this timescale to be studied. As shown in Pap er^ 
significant changes in the PDS of Cyg X-1 do occur on these 
short timescales. The PDS have been rebinned using semi- 
logarithmically spaced bins. 

The normaUza ti on us ed for the PDS is t hat o f 
iBelloni & Hasing eil lll990l) and iMivamoto et"ai] lll992l) . 
where the integral over a frequency range gives the square 
of the fractional root-mean-square (RMS) of that frequency 
interval. In addition, the PDS are plotted in units of frequency 
times power versus frequency (fPf versus /). Together with 
the new data from 2003, the soft state data from Paper l' were 
reanalyzed and corrected for dead-time effects, improving the 
signal-to-noise ratio. Altogether, a total of 614 soft state PDS 
were obtained. 

2.3. Modeling 

iNowakI ilOOd) showed that the PDS of Cyg X-1 in the hard 
state could be successfully fit using Lorentzian compo nents. 
The model was extended by IPottschmidt et alJ (l2003l) . who 
added a power-law component in order to fit the lower frequen- 
cies, especially during the 'failed state transitions'. Building on 
these results, it was shown in Paper I that by introducing a cut- 
off power-law component together with two Lorentzians, we 
were able to provide a good fit to the PDS of Cyg X-1 in all 
states of the source. We could thereby for the first time model 
the complete evolution of the PDS. As the parameters of that 
model differ somewhat from those in the previous studies, we 
present their definition again here. The Lorentzian profiles have 
the form 



Liif) 



HiWiVi 



if - yiY + fWiVi 



(1) 



where Hj is the value of fPf at the peak frequency v,. The di- 
mensionless parameter W, gives a measure of the relative width 
of the profile in the fPj representation. The parameterization 
is chosen such that for a given W,- and Hj, the integral of the 
Lorentzian (and thus its fractional contribution to the total RMS 
variability) is the same for any value of v,-. The cut-off power- 
law necessary in the transitional and soft states is of the form 



Plif)^Af-"e 



(2) 



Where A is the normalization constant, a is the slope, and is 
the turnover frequency. Of the 614 soft state PDS obtained from 
the observations, 143 allowed a fit using all three components 
of the model, one or both Lorentzian components being too 
weak to constrain in the remaining PDS. These fits were first 
performed with all parameters free, and then with the width of 
the first Lorentzian, Wi , frozen to a value of 0.6. While freezing 
Wi did not affect the number of PDS where L2 was detected, 
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it allowed us to better constrain L2 in those cases and the re- 
sults of these latter fits were therefore kept. Table 1 lists the 
results of all fits used in this paper, including the ones from 
IPaper A and is available at the CDS in electronic form only. 
The first six columns give the proposal number and sub-ID 
(e.g., "60090-01-39-00"), the start and stop times in MJD of 
the lightcurve used to calculate the PDS, a number identifying 
the model components present in the fit, and the average 2-9 
keV flux and (9-20 kev)/(2^ keV) hardness during that time. 
This is followed by 1 8 columns giving the parameter value and 
error for the three parameters of the power-law (A, a, fc) and 
the parameters of each of the two Lorentzians L\ and L2 (v,, 
and Wi). 

The model used is not aimed at giving the fullest possible 
description of the PDS of Cyg X-1. For example, we do not 
attempt to model very small features of the PDS, and indeed 
the signal-to-noise ratio is insufficient for us to accurately do 
so. The strength of the model lies rather in its simplicity, al- 
lowing us to accurately fit the main characteristics of the PDS 
in all states, and track changes occuring on short timescales. 
Including a cut-off to the power-law proved essential in allow- 
ing the same model to fit all states, thereby giving a more com- 
plete picture of the power spectral evolution. 

2.4. Definitions of state 

When determining the state of Cyg X-1, it is possible to look at 
characteristics of either radiation spectra or temporal analysis. 
Current definitions of state, derived from either of these char- 
acteristics, are mainly phenomenological. These definitions all 
agree on the 'classical' hard and soft states. However, difficul- 
ties arise when attempting to define the extent of these states, 
and the boundaries of the transitional (or intermediate) state. 
Until the physical processes behind the state transitions are un- 
derstood, these definitions will retain some of their arbitrary 
nature. 

In an extensive study of t he radiation spectra of Cyg X- 1 in 
all states. Iwilms et al!'('2006') showed that there is a continuous 
spectral evolution between the states of the source. While tim- 
ing spect ra are ge nerally able to provide sharper criteria, it was 
shown in'Paper l" that there is a continuum of PDS between the 
hard and soft states, with examples covering all stages of the 
evolution. Therefore, criteria of state such as the temporal fea- 
tures crossing certain frequency boundaries, or the appearance 
of a power-law component, do not by themselves provide good 
indicators of state. However, our previous study also revealed 
that several parameter correlations change behavior at the same 
stage of the evolution, and we therefore regard these changes to 
be a natural marker of the source entering the transitional state. 
Studies have also shown decreased coherence and enhanced 
time lags between the ha rd and soft lightc urves durin g tran- 
sitions, botii in Cyg X-1 dCui et al.lll997bt IPottschmi dt et alJ 
I2OOO') and oflier so urces (such a s GX 339 -4 and XTE J1650- 
500, s ee iNowak et al.. 2002: .Kalemci et al. 2003. respec- 
tively). IgenUoc^^Ol^Qjuse a combination of photon in- 
dex F and time lags between the 2-4 ke V and 8-13keVchan- 
nels to define the state of Cyg X-1. In their long term study of 
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Fig. 1. Relation between the two peak frequencies of the 
Lorentzian components. T he grey points are 1414 hard state 
and transition points from Paper I, while the black points are 
143 soft state points resulting from the analysis in this paper 

the hard state, IPottschmidt et alJ ( l2003l) showed that the state 
is not uniform, and a distinction is made between 'quiet' and 
'flaring' hard states, where flares and 'failed state transitions' 
are more frequent in the latter. The classification of state is thus 
not rigorously defined, and remains an area of active research 
(for a recent review, seelZdz iarski & GierlifiskL2004) . 

For consistency with Pap er J. we have used the same crite- 
ria when defining state, based on the changes seen in the tempo- 
ral components. With this definition, the PDS in the hard state 
only show L\ and L^- As the source enters the transition, the 
power-law component enters the frequency window, and sev- 
eral parameter relations change behavior. In the soft state, is 
not seen or significantly weaker than L\, or the PDS is com- 
pletely dominated by the cut-off power-law component. For 
comparative purposes, we note that our definition is consistent 
with a criterion of the effective spectral inde x (3-12 keV) 
F ^ 2 .4 for the soft state used elsewhere (e.g., IZdziarski et alJ 
l2002h . 

During the 'quiet' hard state, there are periods when a third 
component appears at the higher frequencies of the PDS and 
the hardness increases. We refer to these episodes as the 'hard 
edge of the hard state', and while we d o not mo del them in 
this paper, they are part of the data from lPaoer ll presented in 
Sect. 10 

3. Results 

We now turn to the results of our analysis, starting with the 
study of the soft state data. Following this analysis, we show 
that a decomposition of the PDS is now possible, and present 
the evolution of each temporal component. 

3.). The fyjo Lorentzians 

While the peak frequencies of the two Lorentzian components 
follow a power-law coiTelatio n in the h ard state, the relation in 
the soft state is still uncertain jPaper 1 Fig. 9). As this deeper 
study is partly aimed at resolving this behavior, we begin by 
showing the frequency relation in Fig.^ As can be seen in the 
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Fig. 2. Results of binning in vi (open diamonds) and V2 (solid squares). The lower line is the same as in Fig. 17 of lPaper 1 a fit to 
the hard state points (index 1.20+0.01), and the upper one a power-law with the same index (1.19+0.14) fit to the soft state points 
(binning along V2). Note how the binning along vi deteriorates as the source enters the soft state, but follows the power-law in 
the hard state to the lowest frequencies. The binning along V2 deviates from the hard state fit at the lower frequencies, but extends 
the relation in the soft state to higher frequencies. 



figure, the spread of points is quite large in the soft state (black 
points). To determine whether there is a significant trend also in 
the soft state we will attempt to bin the data points. The prob- 
lem with the binning is the same as in the more common case 
of regression, when choosing which parameter to treat as the 
independent one. In order to determine the effect of the bin- 
ning procedure in this particular case we first performed Monte 
Carlo simulations, where we introduced noise to points follow- 
ing a power-law relation with known index. The resulting dis- 
tributions were then binned along both axes and the indices 
determined. As long as the uncertainties are small, binning the 
data recovers the correct index. However, the results clearly 
show that an increase in noise tends to produce an artificial flat- 
tening effect along the axis of binning. If the amount of noise 
varies between the parameters, the best results are achieved 
when binning is done along the parameter with smallest er- 
rors. Returning to the analogy with regression, this amounts to 
the parameter with the smallest errors being chosen as the in- 
dependent one. In addition, as our frequency range is limited, 
there are border effects as the components approach these lim- 
its, causing a biased distribution in the frequency nearest the 
border 

We now apply the same binning procedure to the data 
points of Fig. [2 In the soft state, the largest spread is along vi. 
Therefore, binning along V2 more accurately reveals the behav- 
ior of the relation. Conversely, for the hard state data, it is the 
spread along V2 that dominates. In this range, binning along vi 
gives a better indication of the behavior of the relation. Figure 
12 shows the result for both binnings. Also included in the fig- 
ure is the power-law fit to the hard state points made in iPaoer L 



with index 1.20 + 0.01, as well as a power-law fit instead to the 
points in the soft state (solid squares). Between the extremes, 
the results of the two binnings are indistinguishable. The flat- 
tening effects seen in Fig.|2]as the components approach the 
frequency limits are in agreement with the results of our simu- 
lations. Our change of binning parameter is thus well founded, 
but the outermost points at either extreme should be treated 
with caution. 

With the better soft state data and the binning along V2, it 
is clear from Fig.|2]that the relation in the soft state follows a 
power-law with an index similar to that of the hard state. The 
result is of course dependent on the where the soft state re- 
gion is assumed to begin. Setting the lower limit too high in 
frequency (v2 <: 10 Hz) results in a large uncertainty as the 
dispersion is high in this region. Lowering the boundary to in- 
clude the asymptotic behavior of the transition (v2 ^ 7 Hz) 
reduces the fitting error, but introduces a systematic effect in- 
creasing the index. For all reasonable choices of this bound- 
ary (i.e., 7 Hz ^ V2 $ 10 Hz) an index of 1.2 is consistent with 
the result of the fit. We will assume a conservative estimate 
of 1.19 + 0.14 as the index for the soft state points, with the 
value derived from the fit with boundary at 8.8 Hz, and the er- 
ror reflecting the variations arising from the range of possible 
boundary choices. The form of the power-law used for these 
fittings is: 

y, = B 1^ (3) 

where Z? is a constant and /3 the power-law index. 

Adopting the fits for the two states, a look at the constants 
of the power- laws reveals the values of B - 0.103 + 0.002 in 
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the hard state and B = 0.19 + 0.05 in the soft state. The ratio of 
the constants is close to 2, and together with the fitted indices 
a shift of the hard state relation by a factor of two is consis- 
tent with the data. Is it therefore possible that the break in the 
relation during the transitions is simply a shift to a higher har- 
monic? In order to accurately model such a shift, three compo- 
nents would be necessary during the transitions. As our model 
only contains two components, it would not catch the change 
if the shift occurs gradually. Unfortunately, attempting to fit 
the PDS with three Lorentzian components shows that a third 
component is impossible to constrain in the data, and rarely 
produces a better fit. We therefore attempt a different approach 
using simulations. 

3.2. Simulations 

In order to test the hypothesis of Li shifting to its first har- 
monic, we conducted a series of Monte Carlo simulations 
where such an effect was included. Artificial power spectra 
were created and fitted with our model of two Lorentzian com- 
ponents plus a cut-off power-law. As input parameter for the 
simulations we use V2, the peak frequency of L2- This is to get 
the same distribution of frequency points as in the observed 
data. For the oth er pararn eters we utilize the parameter cor- 
relations found in 'Paper i We use analytic expressions to ap- 
proximate these behaviors (shown in Fig.O, adding some ran- 
dom Gaussian spread. For each V2 a corresponding vi was de- 
rived assuming the hard state relation between the frequencies 
(lower line in Fig. |2}, and values for H2 and W2 were calcu- 
lated. Finally, a third Lorentzian component, L^, was added 
with v| = 2vi, i.e. the first harmonic of Li. The parameters 
Wi, W^, Hi and were now determined from the relations 
seen in the left panels of Fig. |3] As Wi is measured to cen- 
ter around one value at lower frequencies, and a lower value 
at higher frequencies, these different values were adopted for 
Wi and (solid and dashed lines in panel b of Fig. Irrespec- 
tively). To simulate a shift in power from Li to L| during the 
transitions, the analytical approximation shown in panel a of 
Fig.|3]was multiplied by modulating factors: 



mi(v2) = 



1 



1 + exp[- 



Av J 



m|(v2) = 



1 



l+exp[+^]' 



(4) 



where mi(v2) is the factor for Hi, m|(v2) the factor for H^, 
Vr the value of V2 roughly in the middle of the transition 
(6 Hz $ V, ^ 7 Hz) and Av is a constant related to the width 
of the transition (~ 1 Hz). The values of these parameters are 
constrained by the transition region as seen in Fig.|5] The fig- 
ure shows that the influence of L| must start at V2 ~ 4 Hz and 
completely dominates above V2 ~ 8 Hz. In our simulations, v, 
was set to 6.5 Hz and Av to 1 Hz. 

Finally, a power-law component was added and random 
Gaussian noise was added to the resulting PDS to account for 
the uncertainty of the measured PSD values, and thereby giv- 
ing the same signal-to-noise ratio as in the observational data. 
We use Gaussian noise as the PDS are the results of both av- 
eraging and rebinning, making the fluctuations approximately 
Gaussian. This assumption is least valid at lower frequencies, 
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Fig. 3. The parameter correlations from IPaoer j (gray points) 
and the analytical expressions used to approximate them in the 
simulations (black lines). The two left panels show Hi (panel 
a) and Wi (panel b) as a function of vi. The two right-hand 
panels show the correlations for H2 and W2 versus V2 (panels c 
and d). In panel b, the black line indicates the approximation 
used for Li and the dashed line for L|. Note also the deviation 
from the approximation at the lowest frequencies in panel d. 
These points are from the hard state, when Li approaches the 
lower frequency boundary, and when the hardness is highest. 
As our simulations mainly cover the transition and soft state, 
this deviation does not influence our results. 



as each bin here is the average of only a few points. However, 
these frequencies are not of interest in this context. Examples 
of resulting PDS are shown in the left column of Fig.|4] along 
with the result of the fitting (right column). The difficulty of 
resolving the two components is very clear in the figure. As the 
width of each Lorentzian is very large compared to their sepa- 
ration, they are unresolved even in the pure signal (solid black 
line in the left-hand panels). Thus, even with a large increase in 
signal-to-noise ratio it will not be possible to determine directly 
whether these two components are in fact present in the data. 
This becomes more clear in Fig.|5j which shows the result of 
the fitting routine when Li and L| are equally strong, reflecting 
the most favorable scenario. From the residuals of the fit (bot- 
tom panel), it is clear that a single Lorentzian can successfully 
replace both Li and L|. 

The artificial PDS were put through the same fitting rou- 
tine as the observed data, and the parameters of the result- 
ing fits (right-hand panels in Fig. |4} plotted in the same way. 
Fig. |5] shows the resulting frequency correlation. Fitting the 
artificial PDS with the model used for the data produces re- 
sults nearly identical to those obtained from the observational 
data. Not only is the change in behavior around the transition 
(v2 ~ 4 - 9 Hz) the same, but when comparing with Fig. ^ it 
is apparent that also the general trend of the data points is re- 
produced. The relative dispersion is greater during the soft and 
hard states, with a minimum occuring during the transitions. 
As the components move to lower frequencies, the points in 
Figs, [n and |6lboth show the same gradual increase in relative 
dispersion. However, the dispersion of the soft state points is 
smaller in the simulations than in the data. This difference can 
be traced to the cut-off of the power-law component. As the 
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Frequency [Hz] 

Fig. 4. A sequence of five simulated power density spectra, dis- 
playing the evolution from hard state to the soft state (top to 
bottom). The panels on the left show the simulated model com- 
ponents (dotted lines): a lower frequency Lorentzian Li and its 
first harmonic L|, a high frequency Lorentzian L2 and a cut-off 
power-law. The solid line is the sum of the noise free compo- 
nents, the points indicate the result after the addition of noise. 
The panels on the right show the results of fitting the simulated 
PDS with our model, using only two Lorentzians and a cut-off 
power-law (dashed lines). The solid line is the sum of the com- 
ponents. Note that even when Li and L| are of nearly equal 
strength (e.g., third row from top), they are not resolved in the 
noise free PDS (solid line in left panel), and the sum is well fit 
by one wider Lorentzian component. 

power-law extends through most of the frequency range in the 
soft state, placing the cut-off too high in frequency will capture 
power from L2- Since L2 is weak in the soft state, this artificial 
weakening increases the uncertainty in frequency. In our sim- 
ulations, we can use our knowledge of the power-law to elim- 
inate this effect, thus reducing the dispersion in the soft state. 
When the parameters of the power-law are left completely free 
also in our simulations, the scatter of the soft state points is 
significantly increased. 

It should be noted that the results do not change signifi- 
cantly if the amount of noise introduced is changed, nor does 
it depend to any large extent on the exact expression used to 
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Fig. 5. Result of fitting an artificially created PDS with the 
same model used for the data. The simulated PDS is shown 
in the top panel, with the power-law component as the long- 
dashed line, Li dotted, L| dashed and L2 dot-dot-dashed. 
The middle panel shows the result of the fit using only two 
Lorentzian components, Li (dot-dash) and L2 (dot-dot-dash). 
The bottom panel shows the residuals of the fit. Note that 
there is no trace of the third component. Trying to add a third 
Lorentzian to the fit is not feasible, as the data do not allow it 
to be constrained. 



approximate the parameter relations. Perhaps not surprisingly, 
the results are most sensitive to the modulation function used 
to vary the strengths of Li and L|. This is evident in the H[ - vi 
relation (panel a in Fig. 0- If Li starts to decline while L| is 
still very weak, the fitted Hi values will show a dip at the coiTe- 
sponding frequency, and there is no such feature present in the 
data. We have tested expressions with both asymptotic expo- 
nential behavior and behavior of power-laws of varying index. 
The conclusion is that the rise of L| and the decay of Li must 
be fairly swift. This is also evident from the data in Fig. [l] 
as the transition occurs over a fairly narrow frequency range. 
We found that the observational results were most readily re- 
produced using exponential functions, as in Eq.|4] Using, e.g., 
polynomial modulation functions will achieve good agreement 
with the frequency relation of Fig. ^ but cannot successfully re- 
produce the other parameter relations (gray points in Fig.|3}- In 
particular, the inadequacy in correctly reproducing the parame- 
ter correlations was seen also for modulations of the functional 
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Fig. 6. Results of fitting artificially created PDS with the same 
model used for the data. In the artificial PDS a shift from Li 
to its first harmonic is introduced, successfully reproducing the 
behavior in the data. Note that the general trends of Fig. are 
also reproduced: increased relative spread towards the lowest 
frequencies, low relative dispersion during the transition and 
increased relative spread in the soft state. 



form ( 1 + (v/vo) 1 , typically associated with band-pass fil- 
ters. 

The results from the simulation reproduce the behavior ob- 
served from the data for all parameter correlations. The inter- 
pretation of a shift in power from Li to its first harmonic during 
hard to soft transitions is therefore consistent with the data. 



3.3. Evolution of the PDS 

We now return to the observations, and with the improved re- 
sults from the fits of the soft state, we attempt to separate the 
power spectral components and study their evolution through 
the spectral states. For this analysis we use both the data pre- 
sented here, complemented with the data of IPaoer 1 The goal 
should be to ultimately understand the relation between the 
components of the temporal variability and those of the radia- 
tion spectrum. The transitions are of key importance since they 
are associated with many changes in both timing and radiation 
properties. To measure the relative strength of each PDS com- 
ponent, and its evolution, we computed the fractional RMS by 
integrating over the 0.01 - 50 Hz range. The result is shown in 
Fig- Altogether with the total fractional RMS, a^tot, defined as 



O-tot = 



0"tot 



(5) 



where cr,ot is the total RMS and F the mean flux of the 
lightcurve used in calculating the PDS. We furthermore assume 
that 



(6) 



where cri^y and o-pi correspond to the RMS of both Lorentzians 
and the power-law component respectively. It is assumed that 
both Lorentzians are produced by the same process, indepen- 
dent of that generating the power-law component. Since we ne- 
glect power below 0.01 Hz, the RMS values for the power-law 
component should be seen as lower hmits. 



^liii.** " ■*« ***** II i 
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Fig. 7. Total fractional RMS (triangles) in the 0.01-50 Hz range 
and contribution of each component, plotted against vi . Black 
squares indicate the contribution from Li and L|, grey squares 
from L2 and diamonds from the power-law. Note how the to- 
tal fractional RMS stays nearly constant throughout the hard 
state, dips during the transitions and increases again in the soft 
state. The increase below vi ~ 0.2 Hz coiTesponds to the obser- 
vations with highest hardness ratio, where a third Lorentzian 
component is seen to enter the PDS. These cases make up less 
than 10% of the total hard state observations. 

Figure shows that (Xtot remains quite constant through- 
out the hard state (vi $ 0.65 Hz). The exception to this is the 
range with the highest hardness ratio, i.e. at the lowest frequen- 
cies (vi S 0.2 Hz). In these observations, there is evidence of 
a third Lorentzian component entering the frequ ency window, 
thus increasing a^tot (see, e.g., Nowak,2000: Pottschmidt et alJ 
2003). In the more common hard state, a^tot reaches a plateau 
(0.25 < VI < 0.65 Hz) with a slope of (-0.003 + 0.004) Hz"' 
(zero within Icr, 800 data points). However, the contribution 
from each Lorentzian component is variable. As the power- 
law component enters the window, the Lorentzian components 
weaken. Although the power-law component steadily increases 
in power, it does not fully compensate for the weakening and 
CT^tot declines during the transition. In the soft state the power- 
law starts to dominate the PDS and o^tot once again increases. 

Continuing with the assumptions above, we may consider 
that the components modulate two independent emissions so 
that the flux F can be written as 



F = Fi,, + F, 



pi: 



(7) 



where Fior and Fpi are the fluxes related with respectively 
the Lorentzian and power-law components. However, to derive 
each flux component we need an additional assumption. As 
noted above, the fractional RMS of the combined Lorentzian 
components remains fairly constant during the hard state (^ 
0.28 with a relative dispersion of 4.4%), and shows a decline 
just when the power-law component appears in our frequency 
window in the transitional state. We now assume that 



O'loiV^'lor = const. 



(8) 



for all states, with cr\oi being the RMS of both Lorentzians. 
Thereby associating a smaller fractional RMS to the power-law 
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Fig. 8. Flux relating to the Lorentzian and power-law compo- 
nents as a function of hardness (9-20 keV over 2-4 keV) un- 
der the assumption that the combined fractional RMS of the 
Lorentzians is constant. The Lorentzian components (upper 
panel) follow roughly the same hardness -flux pat tern as the 
source in general (compare, e.g., Fig. 7 of IPaoer j) . However, 
the monotonic decrease of F^i (lower panel) with increasing 
hardness can be explained by the contribution of the power-law 
decreasing. The true hardness-flux behavior associated with the 
power-law is positive, and that is why the spread of points ap- 
pears to systematically decrease. The symbols refer to the dif- 
ferent model components used: two Lorentzians (grey points), 
two Lorentzians and power-law (open squares and black points 
for transitional and soft state points respectively), Lorentzian 
plus power-law (open circles) and only power-law (filled trian- 
gles). 
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Fig. 9. Fractional RMS of the temporal power-law component 
versus flux related with the component. There is considerable 
spread, but the result rules out any strong correlation. The sym- 
bols are the same as in Fig.|Sl 



flux. As the hardne ss-flux relation is known to be positive in 
flie soft state (e.g., IWen et allEoOlt |zdziarski et al. 200^, it 
is likely that the general increase in flux is merely a result of 
the power-law component beginning to dominate the PDS, and 
that the true hardness-flux relation is positive, and the cause of 
the systematic increase in spread. The fractional RMS of the 
power-law (Fig. |9j does not show any strong trend with the 
associated flux, and the results are consistent with a constant 
ratio. As seen in the figure, the dispersion is greatest when the 
power-law component is weak, and decreases as the compo- 
nent becomes stronger We have calculated the median values 
in intervals, and they all coincide (within 2cr) with the median 
when the power-law is the only component seen in the PDS 
(black triangles in Fig.|5). This value is 0.24 ± 0.004, and the 
sample standard deviation is 0.03 (14%). 

4. Discussion 

In order to establish whether a shift of Li to L| is truly a possi- 
ble scenario, all parameters must be examined. We also discuss 
the results of the PDS decomposition, which provide clues to 
the relation between the components. 



component, we can determine the individual flux contributions 
in Eq.0 These quantities are shown plotted versus hardness in 
Fig- 13 and the RMS of the power-law (o-pi) is plotted versus 
flux related with that component in Fig.|5] 

As can be seen by comparing Fig. |8] with the total flux- 
hardness relation in C yg X-1 (e.g.. Fig. 7 in Paper 1; Fig. 4 in 
IZdziarski et alJ2003) . it is the behavior of the Lorentzian com- 
ponents that appears to determine the flux-hardness correlation. 
Furthermore, the flux related with the Lorentzians is reduced 
in the soft state. The weakening of the Lorentzians in the soft 
state is thus due to an actual weakening of the corresponding 
emission component, and not simply due to another emission 
component growing stronger. The flux related to the power-law 
component decreases monotonically with increasing hardness. 
Note however that the spread also increases systematically with 



4.1. Shift to first harmonic 

If the apparent change in the frequency coiTelation is caused by 
L\ shifting to a higher harmonic, this should be visible in the 
other parameter coiTelations as well. The behavior of H\, W\ 
and hardness ratio as a function of fitted Lorentzian frequency 
does indeed show some clear changes during the transition. 

As noted above, the value of the width parameter W\ ap- 
pears to center around 1.2 in the hard state and another, lower 
value of ~ 0.7 in the soft state (panel b in Fig. I^l. While this 
can simply indicate that L\ becomes more defined in the soft 
state, it is interesting to note that a shift to twice the frequency 
also predicts an increase in the coherence of the signal by a 
factor of two. While W\ is not directly a measure of the co- 
herence, the decrease is consistent with the higher coherence 
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expected if observing Li in the hard state and L| in the soft 
state. Our parametrization does not directly give a value for the 
more generally quoted quality factor Q - v/Av. The relation 
between W and Q is: 



W = 2 



1-1 + 



1 

4fi2 



(9) 



Another clue may be obtained when looking at the behav- 
ior of H[ (panel a in Fig. |3}- There is considerable spread in 
the hard state, but as the source approaches the transition, the 
value appears to increase slightly, only to decrease sharply as 
the transition is reached. In the context of a shift to the first har- 
monic, the rise would start when the strength of L| becomes 
comparable to that of Li, and the rapid decline is seen when 
Li begins to weaken. As noted in Sect. 13.21 above, this has to 
happen fairly rapidly, consistent with the sharp drop in Hi . 

As the hardness is not a parameter of either Lorentzian 
component, it gives an independent measure of the evolution 
of the source. It also allows a comparison to determine if the 
shift in the frequency relation is related to changes in vi or 
V2 (or both). Plotting the hardness against these two parame- 
ters ("Fig. I lot does reveal some differences. For V2, the relation 
is a smooth power-law, with a flattening at the highest hard- 
ness. However, while the vi relation also describes a power-law, 
there is what appears to be a shift during the transition. Though 
not immediately evident, this shift indicates that the change in 
the frequency relation is indeed caused by changes in vi, and 
is furthermore what would be expected if a shift to the second 
harmonic occurs. 

Taken together, the above features provide rather convinc- 
ing evidence for the idea of a shift from Li to L| during the 
transition. The observed properties are certainly consistent with 
that interpretation. Any suggested model needs to explain all 
the features consistently, and we find that invoking the shift to 
L| is a simple but successful scenario. 

We al so note that in their long -term study of the PDS of 
Cyg X- 1 . IPottschmidt et alJ ( l2003b report a shift in power be- 
tween their Lorentzian components. In the hard state, the PDS 
are dominated by Li and L2, while during the transitions and 
'failed state transitions' Li decreases in power and L3 strength- 
ens. As their identification of components is not the same as 
that used here, a direct comparison is not possible. However, 

the mechanism sugges ted above is consistent wi th their results. 

' 1 I — ^ — • 

In addition, results by Pottschmidt et al. (2005) from interme- 
diate state PSD indicate that Li may decrease earlier at higher 
energies. 

An appearance of a QPO with harmonic content has been 
reported by Belloni et al. (2005J in GX 339-4 in connection 
with the source shifting from the hard state to the intermediate 
state. If the change in the frequency relation observed in Cyg X- 
1 during state transitions is indeed caused by L| growing in 
power, the question arises if the appearance of a harmonic is a 
more general feature of the intermediate state of black hole bi- 
naries. The appearance of (or in the case of Cyg X-1, a gradual 
shift to) a harmonic could then be used as an indicator of spec- 
tral state. Certainly the change in the frequency relation may 
used as a gauge of the state transition in Cyg X-1, and may 
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Fig. 10. Hardness ratio as a function of vi and V2 respectively. 
While the relation appears fairly smooth for both components, 
there is evidence of a shift in the vi relation. The filled squares 
show the results when binning the points along frequency. The 
dashed lines show the behavior of the v,-HR relation. For V2 
(lower panel), a single power-law describes the relation through 
all the states. However, for vi it is necessary to shift the power- 
law during the transition, as seen in the upper panel. 



as such be added to the growing list of possible ways to deter- 
mine the state of the source. It is likely that all these indicators 
are responding to the same changes in the source (e.g., accre- 
tion geometry) rather than being directly related. Nonetheless, 
when combined they do provide constraints on proposed mod- 
els for the state transitions. 

In the cont ext of the relativistic p re cessio n model 
discussed in IPaperj JStella & Vietril Il998l Il999t 
IStella. Vietri.&Morsinklll999h V, ^ 2v„" and V2 = v, 



peri 

the 



where Vnod is the frequency of nodal precession and Vper 
frequency of periastron precession. The shift would then mean 
that v| - 4vnod. Assuming the frequencies are produced close 
to the inner edge of an accretion disk (or picked out, e.g ., 
from a narrow transition region, see IPsaltis & Normanl200rit) . 
extending the relation to higher frequencies allows tracking of 
the disk to smaller radii. Assuming a prograde scenario with 
dimensionless specific angular momentum a, = 0.49, 



R„ 18 Ma/ \10Hz; 



llOHz/ ' 



(10) 



10 



M. Axelsson et al.: Probing the temporal variability of Cygnus X-1 into the soft state 



where Ri^ is the inner disk radius, Rg the gravitational radius, 
and M the black hole mass. The dependence on |a,| is very 
weak in the studied range. 

4.2. Decomposed power spectrum 

Although the RMS of the two Lorentzian components show 
different behaviors and evolution during the hard state, the sum 
remains constant. While this is no direct evidence that the com- 
ponents are tied, it is improbable that two unrelated compo- 
nents would combine to a constant total. We therefore view the 
behaviors of the Lorentzians as indicative of them arising from 
the same physical process. One possibility (especially in the 
context of the relativistic precession model discussed above) is 
a scenario where the two components share the same, constant 
energy. As the amplitude of one grows, the energy available to 
the other will drop and vice versa. 

An extension of this reasoning indicates that the variabil- 
ity giving rise to the power-law component is instead indepen- 
dent of the process behind the Lorentzians. When the power- 
law enters the frequency window, the total RMS drops, reach- 
ing a minimum when the PDS c omponent s are a ll roughly of 
the same strength. As shown bv IZdziarskil ( l2005h . this is pre- 
cisely what is expected if the measured flux arises in different 
and unrelated processes. We therefore conclude that there are at 
least two unrelated radiation components producing the PDS of 
Cyg X-1. One dominates the hard state and provides the vari- 
ability seen in the Lorentzian components, the other dominates 
in the soft state and produces the power-law component. The 
dr op in total RMS during tr ansitions can also be seen in Fig. 3 
of lPottschmidt et alJ ll2003l) . 

This scenario fits with previous analyses of the radia- 
tion spectrum, which typically include both a soft black- 
body component and a harder compon ent usually ascribed to 
Comptonization in a hot corona (e.g., I^erliiiski et al."l997', 
[oierliiiski et al. 1999). Recent st udies of the energy-dei^ endent 
variability of Comptonization dGierliiiski & Zdziarskil I2OO5I) 
show promising results in tying together temporal and spectral 
components. Work along these lines has also been presented by 
iMivamoto et all (Il994i) for the sour ce G S 1 124-683, and more 
recent ly bv IVienarca et alJ ( l2003l) and IZvcki & Sobolewskal 
il200.4 . 

An interes ting possibility for the t wo components has been 
suggested by Ibragimovetal ] ^005). They find that the ra- 
diation spectrum of the hard state of Cyg X-1 can be well 
fit using both a thermal and non-thermal distribution of elec- 
trons, associated with an inner hot corona and magnetic flares 
above the accretion disk respectively. The transitions into the 
soft state can then be explained by the inner edge of the disk 
moving to smaller radii. As this happens, the corona is reduced 
(thus weakening the thermal Comptonization component), and 
emission from the non-thermal distribution connected to the 
flares becomes more prominent in the radiation spectrum. It is 
tempting to associate the Lorentzian components with a transi- 
tion region between the cool disk and the inner quasi- spherical 
corona, and similarly tying the power-law component of the 
PDS to emission from a non-thermal distribution of electrons, 



associated with the flares above the disk. Not only does the 
relative strength of the power-law increase, but also the cut-off 
moves to higher frequencies as the source enters the transitional 
and soft states, when the disk is expected to extend closer to the 
compact object. Comparing the results of this study with those 
gained fro m fitting the r adiation spectrum of Cyg X-1 in all 
states (e.g. Jlbragimov et al.i.2005i: .Wilms et al..2006.) indicates 
that there is no trivial, direct connection between the compo- 
nents of the radiation spectrum and those of the PDS. In order 
to build a complete physical picture, the behavior of the radio 
jet of the source also needs to be included. Studies have found 
correlations bet ween the radio and X-ray flux in Cyg X-1 (e.g., 
iBrocksopp et al.>.1999: .Poolev et al...l999.), and the radio flu x 
has been con nected to the state tra nsitions dWilms et alJ2006l) . 
Furthermore, Migliari et al. (2005) have found correlations be- 
tween radio luminosity and temporal features of the PDS in a 
number of X-ray binaries. 

If the power-law PDS component is taken to have constant 
fractional RMS, the value cannot be the same as that for the 
Lorentzian components. While the distribution in Fig. |9l does 
not rule out a trend in the fractional RMS, drastic changes 
are not consistent with the data. We have performed Spearman 
rank correlation analysis and find only very weak coiTelations 
between the parameters of the power-law and those of the 
Lorentzians. If the power-law component is connected with 
simple shot noise, the flat behavior seen in Fig. |9l indicates 
that the increase in flux is due to the amplitude of the shots 
increasing, and not to a greater number of shots, as that would 
reduce the observed fractional RMS. The results are consis- 
tent with th e RMS being proportional to the flux, as shown by 
lUttlev fc McHardv t200li) . 

In a study of the relation between RMS and flux on both 
longer and shorter timescales, iGleissner et alJ ( l2004l) find that 
the fractional RMS in the 1 - 32 Hz frequency band to be con- 
stant throughout the hard state changes of PDS, similar to the 
results presented above. They investigate the RMS-flux corre- 
lation in terms of a linear relation and show that such a relation 
is present in both the hard and soft states of Cyg X- 1 . However, 
the slope and intersect of the linear relation have different val- 
ues in different states. While the paper mainly concentrates on 
interpretations in terms of a component with either constant 
flux or constant RMS, the authors also mention the possibil- 
ity of a component with constant fractional RMS. We find this 
interpretation to fit well with the results presented in this pa- 
per, where both the sum of the Lorentzian components as well 
as the power-law component display such behavior. As noted 
above, the two cannot have the same value for a constant frac- 
tional RMS, and it is therefore natural that the parameters the 
RMS-flux relation be different in the hard and soft state, as the 
power-law is not visible in the former but dominates the latter. 
During the intermediate state, both components are present, re- 
sulting in parameter values in-between the h ard and soft state 
values, as reported bv lGleissner et al.lll2004l) . 

5. Conclusions 

We have conducted an in-depth investigation of the soft and 
transitional state PDS of Cyg X-1 both reanalyzing previous 
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data and adding data from new soft state observations. Fitting 
the PDS with a model of two Lorentzian components and cut- 
off power-law, we are able to extend the frequency relation 
dWiinands & van derKhsll999tlPaDer]h into the soft state. We 
show that the index of the relation is consistent with that found 
in the hard state. This further supports the identification of the 
frequencies as the relativistic precessional frequencies. 

Using an approach of Monte Carlo simulations we show 
that a shift in Li from vi to its first harmonic v| = 2vi is con- 
sistent with the behavior of all parameters. We conclude that 
this scenario provides a simple explanation of the observational 
data. 

Combining these results with the hard and transitional state 
results from Paper 1, we decomposed the PDS and studied the 
evolution of the individual components. The results show that 
although each Lorentzian is variable, the total fractional RMS 
is constant throughout the hard state of the source. In contrast, 
when the power-law component enters the frequency window, 
the total RMS drops. Together with the lack of correlation be- 
tween power-law and Lorentzian parameters, this suggests that 
these components arise in different regions and/or through dif- 
ferent processes. Using the observed behavior, we attempted to 
derive the relation between hardness and flux for the PDS com- 
ponents, as well as the relation between flux and RMS for the 
power-law. We conclude that the weakening of the Lorentzian 
components as the source enters the soft state is due to an actual 
decrease in the associated flux component. 

The results show that at least two variability components 
must be present in the 2-9 keV energy range studied here. In 
order to better understand both Cyg X-1 and similar objects, 
studies of both temporal and radiation spectra need to be taken 
into account. Only by combining the insights from these fields 
can we hope to understand these sources. 
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